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Recent direct H(z) data indicate that the parameter H{z) may wiggle with re- 

bc , 

, spect to z. On the other hand the luminosity distance data of supernovae flatten the 

o 

CN ; wiggles of H{z) because of integration effect. It is expected that the fitting results 
>v 

d ■ can be very different in a model permitting a wiggling H(z) because the data of su- 



pernovae is highly degenerated to such a model. As an example the natural phantom 
dark energy is investigated in this paper. The dynamical property of this model is 



studied. The model is fitted by the direct H(z) data set and the SNLS data set, 
respectively. And the results are quite different, as expected. The quantum stability 
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of this model is also shortly discussed. We find it is a viable model if we treat it as 
an effective theory truncated by an upperbound. 



I. INTRODUCTION 

the universe is one of the most significant cosmological discoveries 
J]. The decisive evidence of the present acceleration is witnessed by 
supernovae. The principle of this conclusion is based on the fitting of LCDM (cold dark 
matter with a cosmological constant) model by the data of the luminosity distance of the 
supernovae. The phenomena of this cosmological acceleration is of very interest, for which a 
large number of models have been proposed besides LCDM, and several of them have been 



fitted by luminosity distances of supernovae, for a review, see J2]. 

In a cosmological model, one calculates the luminosity distances as follow, 

n l + z [' H , 

where H(z') denotes the Hubble parameter, and Ho represents its present value. Then one 
can constrain the parameters in the model by using the observation data of supernovae 
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z 


0.09 0.17 0.27 0.40 0.88 1.30 1.43 1.53 1.75 


H(z) (km s^Mpc" 1 ) 


69 83 70 87 117 168 177 140 202 


68.3% confidence interval 


±12 ±8.3 ±14 ±17.4 ±23.4 ±13.4 ±14.2 ±14 ±40.4 



TABLE I: The direct observation data of H(z) 



6J] (see [7j also). 



through x 2 ° r other method. A deficiency of this method is that one obtains luminosity 
distance through an integrating to Hubble parameter H, and therefore, the fine structures, 
such as wiggles on H, can not show themselves in such a method. For example, compared 
with 

H{z)/H = 1, (2) 

H(z)/H = 1 (3) 
1 + sm{nz) 

is surely a different model. But for a large n, they always share the same confidence region 
in fittings by using luminosity distances data. And different n (for large n) also share the 
same confidence region. We see that some types of fine structures of H is highly degenerate 
to the luminosity distance data. To break this degeneration one needs the observational 
data of H(z), not only an integration of H~ l (z). 

Fortunately there is a newly developed scheme to obtain the Hubble parameter directly 
at different redshift [3( , which is based on a method to estimate the differential ages of the 
oldest galaxies [4]]. By using of the previously released data 5], Simon et al. obtained a 
sample of direct H(z) data in the interval z G (0, 1.8) p, almost as the same interval of the 
data of luminosity distances from supernovae. We show this sample in table I. 

Table I displays an unexpected feature of H(z): It decreases with respect to the redshift 
z at redshift z ~ 0.3 and z ~ 1.5, which means that the total fluid in the universe behave 
as phantom. This information of dynamical property of the universe is very difficult to be 
drawn from the data of supernovae. This feature of H(z) implicates that the dark energy 
component of the cosmic fluid must behave as phantom sometime, which can be proved by 
the following argument. In standard general relativity for a spatially flat universe, which is 
implied either by theoretical side (inflation in the early universe) ,or observation side (CMB 
fluctuations [8|]), the Friedmann equation reads, 

H 2 = ^~Ap m + Pde), (4) 



3fi 
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where p m denotes the density of dust matter, p de stands for the density of dark energy, and 
p represents the reduced Planck mass. Differentiate with respect to the redshift z, we derive 
from (Hj), 

2H— = 1 ( dPm I dpde ) (5) 
dz 3/i 2 dz dz ' 

where a prime denotes derivation with respect to z. Clearly, if dH /dz < at some redshift 
(as shown in table I), one concludes dpd e /dz < since dp m /dz > 0, which means the dark 
energy behaves as phantom. 

The present (or at very low redshift) phantom behavior of dark energy is also implied by 
the supernovae data [9]. Generally speaking, a simple phantom field (scalar field with kinetic 
term of false sign) is quantum mechanically unstable. However, several evidences imply that 
our present 4 dimensional standard model and general relativity is not the final theory. The 
phantom model can be treated as reduced theory of more fundamental theory, in which there 
is no field behaves as phantom [2]. Thus the stability problem may be evaded. Actually, 
many of the reduced theories do contain phantoms, as the ones coming from string and/or 
M-theory compactification, or higher-derivative supergravities, or modifications of Einstein 
gravity itself, for example, such a field may be motivated from S-brane constructions in 
string theory [h]]. Moreover, there exist examples in which an effective phantom and/or 
quintessence description of the late time universe naturally emerges, even when the starting 
theory does not clearly show the phantom and/or quintessence structure ll|. Therefore 
it may be reasonable to investigate such models as an effective theory. Phenomenologi- 
cally, the cosmological models with phantom matter have been investigated extensively [12 ]. 
Also, urged by observations, the models with dark energy whose EOS crosses —1 have been 
investigated in [3]. 

However, the H(z) data in table I implies the EOS of total fluid in the universe crosses 
— 1, not only the dark energy sector. Moreover, the phase oscillation over deceleration phase 
and acceleration phase is clear through the history of the universe. In a fitting in frame of 
LCDM model, the point z ~ 1.5 is determinately beyond 1-a level jg]. Furthermore, it is 
shown that the data point near z ~ 1.5, which dips so sharply and stays clearly outside of 
the best-fit of the LCDM, XCDM and 0CDM models studied in fl. Contrarily, a study 
show that the model whose Hubble parameter is directly endowed with oscillating ansatz by 
parameterizations fit the data much better than those of LCDM, IntLCDM, XCDM, IntX- 



CDM, VecDE, IntVecDE 



14J. However no previous physical dark energy models possessing 
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this oscillating property. Therefore, it deserves to present a physical model in which the 
EOS of total fluid crosses — 1. 

In this paper we put forward a model in which phantom field with natural potential ,ie, 
the potential of a pseudo Nambu-Goldstone Boson (PNGB), drives the universe. We shall 
show that in such a model all the features of H in table I can be realized naturally, and the 
fitting results of the parameters in this model are rather different according to supernovae 
and direct H(z) data. PNGB is an important idea in particle physics. It emerges whenever 
a global symmetry is spontaneously broken. There are two key scales of PNGB generation. 
One is the scale at which the global symmetry breaks, denoted by /, and the other is the 
scale at which the soft explicit symmetry breaks, denoted by C. Under this assumption the 
potential of PNGB reads, 

V = C 4 ( 1 ± cos(^)Y (6) 



f ' n 

Inflation model driven by a scalar with such a potential was firstly studied in |15fl. Generally 
speaking in the context of inflation model the cosine function in potential never completes 
a cycle. The scalar PNGB can also play the role of dark energy 16| . In this scenario 
the energy scale of the the global symmetry breaking / keeps about the same as the case 
of inflation model ,ie, the Planck scale. Contrarily, the scale of explicit symmetry breaking 
decreases to an extremely low scale, ie, 10 _3 eV, which is comparable to neutrino mass yielded 
by Mikheyev-Smirnov-Wolfenstein (MSW) mechanism. Phenomenologically, the natural 
potential has been generated to solve the coincidence problem, in which the cosine function 
in potential oscillates many cycles 17], and therefore the densities of dark energy and dust 
can be comparable several times in the history of the universe. But the previous models with 
PNGB dark energy can not realize the feature that the EOS of total fluid in the universe 
crossing — 1. This feature appears naturally in the present phantom natural dark energy 
model. 

In the next section we shall present the phantom natural dark energy model and inves- 
tigate some dynamical properties of it. In section III, we fit this model by using the SNLS 
data and direct H(z) data, respectively. The main conclusions and some discussions appears 
in the last section. 
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II. THE MODEL 

We work in the frame of the standard 4 dimensional general relativity. The phantom is 
characterized by a false sign of kinetic term in the Lagrangian, 

£,= ^-V(^ (7) 

where and in the follow, we take the signature (— , +, +, +). In the present model a phantom 
field with generalized natural potential plays the role of dark energy. In an FRW universe, 
Pde in PJ becomes 

P^—^ + V^), (8) 

and the pressure of the scalar reads 

P* = ~^-V(0), (9) 
where a dot denotes derivative with respect to time. And the equation of motion of <fi reads, 

-0 + 3#0+^ = O. (10) 
dip 

Based on the former researches, we phenomenologically generalize the natural potential to 
the following form, 

V((f>) = V (l + Acos(p£)). (11) 
With the new dimensionless variables below, 

* 4 vfe (12 > 



" 4 vfe- (13) 

' 4 7& (14) 

A Wo 



b 4 (15) 
the dynamics of the universe can be described by the following dynamical system, 

x' = ~x{-2x 2 + I 2 ) -3x + ^.py/AW - (y 2 - b 2 ) 2 , (16) 

y' = \y{-2x 2 + I 2 ) - y^ pX y-^AW-(y 2 -b 2 y, (17) 

I' = i(-2x 2 + / 2 )-i, (18) 

b' = ^ b (-2x 2 + l 2 ) } (19) 
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where a prime stands for derivation with respect to s = ln(l + z). Note that the 4 equations 
(flBjl . ( JT71) . ( fl8j) . ( fl~9l) of this system are not independent. By using the Friedmann constraint, 
which can be derived from the Friedmann equation, 

-x 2 + y 2 + l 2 = 1, (20) 

the number of the independent equations can be reduced to 3. There are four critical points 
of this system satisfying x' — y' = V = b' = appearing at 



' / o. n i. h ±\/y^; (21) 



, / 0. yy: I. h ±^I-L_. (22) 

All of them satisfy the Friedmann constraint (1201 . To obtain real values of the variables 
at the singularities we see that if A > 1 only the former two exist, if A < —1 only the 
latter two exist, and only for —1 < A < 1 all of the four critical points exist. The critical 
points imply that the universe will enter a pure dark energy phase at last, if the singularity 
is stationary. To investigate the properties of the dynamical system in the neighbourhood 
of the singularities, impose a perturbation to the critical points, 

5x' = En6x + E 12 5y + E u 5b, (23) 

by' = E 22 5y + E 24 5b, (24) 

51' = E 33 6l, (25) 

5b' = 0, (26) 

where we have used (l2Tj) or ff22l . and the components of the eigenmatrix reads, 

E n = -3, (27) 

E 12 = -y/6py( y 2 -b 2 )a- 1 , (28) 

E u = VQp[A 2 b 3 + (y 2 - b 2 )b]a~\ (29) 

E 22 = VQxp{y 2 -b 2 )a-\ (30) 

E 2A = -VQpxy^^b 3 + b{y 2 - & 2 )]a _1 , (31) 

E 33 = -3/2, (32) 



where 



a = ^/A 2 b 4 - (y 2 - b 2 ) 2 . (33) 
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The 4 eigenvalues of this linear system reads 

\ 1 = -3, A 2 = V6xp(y 2 - b 2 )a'\ A 3 = -3/2, A 4 = 0. (34) 
The property of A 2 is rather complicate around the singularities. For example, it goes to 



different values along different pathes around the singularity x — 0, z = 0, b = 1 /yl + A . 
The 6 repeated limits read, 

(35) 



hm hm hm A 2 

^n.r^n 6 , i 


= o, 


lim lim lim A 2 

:T ^n t , i 


= 


lim lim lim A 2 

2->nf, . i .r-+n 

" Vl+A 


= o, 


lim lim lim A 2 

, , i ^^n.T-^n 

° Vl+A 


= o, 


lim lim lim A 2 


= o, 


lim lim lim A 2 


= ip 



3A 
1 + A 



3A 
l + A 



(36) 
(37) 
(38) 
(39) 

(40) 



Hence the limit of A 2 does not exist at the singularities. However, we see that the real 
parts of the limits keep zero independent of pathes, which means that the system reaches an 
indifferent equilibrium. In such a de Sitter universe at the critical point the kinetic energy 
of the phantom and dust matter vanish, but the potential energy can reside at any values, 
which depends on the initial values of kinetic energy, potential energy, dust density and the 
Hubble parameter. 

As we have pointed out in section I, the data of supernovae is insensitive to the oscillating 
behaviour of H(z). In this section we show the fitting results by the direct H(z) data and 
SNLS data by x 2 -statistics, respectively. The H(z) data have been used to constraint models 
in [3]. Here we adopt SNLS data [l^ . which is believed to be more consistent with 

CMB data. Figure 2 displays the fitting results. We set A = 1, which means we adopt the 
original PNGB potential, b(z = 0) = 0.616, <f)/fi(z = 0) = 0.022, H Q = 72km s" 1 Mpc" 1 



20| | . In figure 2 we find an extraordinary property of (a): the 68.3% confidence contour is 
disconnect. The physical explanation is that the data set of direct H(z) is too small, that 
is, the data do not distinctly illuminate how many "wiggles" inhabit on H(z). New wiggles 
may hide in the gaps of the data set, which leads that a much bigger p lies in the same 
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20 
18 
16 
p 14 
12 
10 




0.28 0.3 0.32 0.34 0.36 0.38 0.4 



0.28 0.3 0.32 0.34 0.36 0.38 0.4 



FIG. 1: The fitting result of the parameters VL m (/q) and p. (a) The 68.3% confidence contour 
plot by using the direct H[z) data in table I. (b) The 68.3% confidence contour plot by using the 
SNLS data. 




FIG. 2: The deceleration parameter q with best fit Q m = 0.397 and p = 11.0. 

confidence region as a smaller p. (b) clearly shows that the resolution of supernavae data is 
very inefficiency to the oscillating behaviour of H(z). We show the deceleration parameter 
q in figure 2 with best fit values of Q m and p by direct H(z) data. Figure 2 illuminates that 
the universe oscillates between deceleration phase and acceleration phase. 
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III. QUANTUM STABILITY 

A severe problem of any phantom field is quantum stability. In practice, we do not require 
that the phantom is fundamentally stable, but quasi-stable, which means, its lifetime is larger 



than the age of the universe. This problem has been discussed in 21[ 22]. Here we follow 



the investigations in 2JJ . The simplest interaction between phantom and graviton takes the 
form, 

-> h + 01 + 02 , (41) 

where h denotes the gravitational fluctuations on an FRW background, 0i, 02 represent other 
two phantom fields. We note here that, though h will be different if we take a Minkowski 
background, but, the difference is tiny and negligible in the spacetime region we considered 
for this interaction. Here we consider a series expansion around the initial value of the 



numerical example we studies above, <\>{z = 0)//i = 0.022. Based on the discussion in 21], 
we set the interaction term, 



C t = I(^)ly'"(0 o ) 

jj 31 



i3 



A\ e (ixh)<j) 3 , (42) 



Where A P is defined as 



Ae = R—i sm G°— ) ~ 10 > ( 43 ) 

where p takes the best fit value in last section p = 15.3. Clearly, if we treat phantom as 
fundamental theory, it will be unstable, and the reaction rate goes to infinity because the 
volume of the phase space of 0, 0i, 02 goes to infinity. However, if we treat it as an effective 
theory which is only valid below some energy scale A, the reaction rate F becomes, 

r ~ \l— , (44) 

where the effective mass of the phantom is defined as 

m, ^ (-AV") 1 / 2 ~ pfiJ^ = 10- 6 V4 1/2 - (45) 

Here and the following, we take A = 1 without special announcement. The phantom field as 
an effective field is viable if its reaction rate T is smaller than the present Hubble parameter, 

r < H ~ 10~ 6 V, (46) 
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which means, 

A 2 — < Ho, (47) 

that is A < 10 58 /x. In fact, any effective theory is valid at such a high energy scale, which 
is far beyond our present lab energy scale, surely be a perfect effective theory. But, besides 
the decay channel as shown in (14ip , we must consider the cases that one phantom decay into 
sever al p articles. By summing over all these possibilities, one arrives at the total reaction 



rate 



2l| 

A 



tot 



r 



1 - 



21 - 2 



(48) 



We see that if A is smaller than fi, T tot will keep the same order of T, which is not a stringent 
constraint. 

However, as an effective theory, one must include all possible terms compatible with the 
symmetry of the Lagrangian up to finite orders to guarantee the renormalizablity of the 
theory-contributions from high order terms are much suppressed which we can neglect up 
to required precision. A most famous effective theory is four-fermion interaction theory, as 
the low effective theory of electro-weak interaction. Expands the propagator of the gauge 
boson in electro- weak according to the mass of W boson, we see that the derivative coupling 
appears. Hence a derivative coupling in an effective theory is quite reasonable. We consider 
the interaction Lagrangian of graviton and phantom, with an approximate global symmetry, 

^ = -^M(d^d4>)? , (49) 

where 7 is a constant of order 1. The reaction rate becomes 

7 2 A 6 



(50) 



which should be smaller than the present Hubble parameter. Therefore we reach 

A 6 < H ° m f* ~ 1(T 12 V 6 . (51) 

7 

The key difference between our result and the result in [2l|] dwells at the effective mass of 
the phantom field. In fact, our result the reaction rate T in (loTJj) is is suppressed by a factor 
p. On the observational side, we see that the 1-a confidence region form a confidence tower, 
no clear upper bound of p. On the theoretical side, one hardly find principles to constrain 
p in the natural potential in this cosmological context. Physically, the effective mass of the 
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phantom in the present model can be notably larger than H = 10 eV, which is taken as 



the mass of the phantom in 211 ] . For example, if = 10 -27 eV, which is quite beyond our 
present abilities of accelerators, A will exceed ITeV and hence it is also beyond our present 
lab energy scale. Thus, the present model is promising due to this suppress mechanism for 
derivative coupling. Also, we note here that a larger A is helpful to increase the mass of the 
phantom, which can be seen from ( |45|) . 

IV. CONCLUSIONS AND DISCUSSIONS 

To summarize, this paper illuminates that direct H(z) data is much more efficient than 
the supernovae for the fine structures of Hubble diagram. 

We first put forward a model based on the previous studies on the PNGB. In this model 
the total fluid in the universe may evolve as phantom in some stages, which contents the 
direct H(z) data in table I. Then we study its dynamical properties and find its critical 
points. And we also study the stability about the singularities of this system. 

In section II we fit our model by using H (z) data and supernovae data, respectively. The 
results are quite different, as we expected. Because the sample of H (z) data is too small, the 
confidence contour is disconnect, which means that we still lack enough information about 
the oscillations of H(z). We hope the future observations offering more data of H(z) so that 
we can investigates the history of the universe in a more detail way. 

In section III we investigate the stability of the present model. Our treatise is to treat 
the phantom model as an effective model truncated at some energy scale A. As the previous 
studies, we find that the coupling to graviton needs a truncate scale much larger than the 
lab energy scale, if we require the lifetime of the phantom is longer than the universe. 
Different from the previous studies, we find that the derivative coupling between phantom 
and graviton is viable due to the special potential of the present model. 
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